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I. INTRODUCTION 

Since the advent of the renormalization group (RG) 
techniques pj-lj], a lot of effort has been put into the 
calculation of the /?-functions and anomalous dimensions 
of specific theories. For instance, full-fledged two-loop 
formulae for non-supersymmetric gauge models became 
available as early as in 1984 thanks to the seminal viforks 
by Machacek and Vaughn [5|-l7| . In the context of super- 
symmetry (SUSY), the need to adopt the existing ma- 
chinery for the soft SUSY-breaking sector postponed the 
arrival of the first generic two-loop results for about ten 
years J8l-ll5l|. Since then, there have even been attempts 
to go beyond two loops in the literature, c.f., |16l.ll7l|. 

For the sake of simplicity, in many of the pioneering 
works the gauge group was assumed to contain at most 
one abelian gauge factor. The point is that with more 
than a single gauged U{1) in play, a new qualitative fea- 
ture requiring a dedicated treatment emerges. This is 
due to the fact that abelian field tensors F^y are not 
only gauge-covariant but rather gauge-invariant quanti- 
ties and, thus, unlike the non-abelian ones, they can con- 
tract among each other without violating gauge invari- 
ance, giving rise to off-diagonal kinetic terms 



the tree level Lagrangian at a certain scale, they are in 
general re-introduced by the renormalization-group evo- 
lution b 
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Moreover, even if such terms happen to be absent from 
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The reason is that the anomalous di- 
mension 7 driving the relevant renormalization group 
equations (RGEs) are in general non-diagonal symmet- 
ric matrices in the gauge-field space, thus giving rise to 
off-diagonal corrections to the gauge boson propagators. 
Tliese, in turn, require extra counterterms in order to 
retain renormalizability. 

Actually, there are several exceptions to this basic rule. 
For instance, it can be that all the relevant U{1) couplings 
originate from a common gauge factor and, thus, barring 
threshold effects, all of them happen to be equal at a 
certain scale. In such a case, accidentally, the charges 
and the gauge fields can be simultaneously rotated at the 
one-loop level so that no off-diagonalities pop up in 7 
[lot |20| and one can use the simple form of the RGEs 
for individual gauge couplings. This is relatively easy to 
implement in the non-SUSY case where only the gauge 
sector has to be taken into account; the only price to be 
paid is the presence of continuous charges in the game. 

In supersymmetry, the U{1) gaugino soft masses can 
also mix, and thus one has to deal with the non- 
diagonalities in the gaugino sector too. Again, the ro- 
tated basis can be helpful if both , gaugino masses and 
gauge couplings, unify at the same scale. However, this 
method is consistent only at the one-loop level where the 
evolution equations for the gauge couplings and gaugino 
soft masses essentially coincide. At two loops, Yukawa 
couplings and trilinear soft SUSY breaking couplings en- 
ter and the relevant algebraic structures are independent 



of each other which, in turn, renders this approach use- 
less. 

For the non-SUSY gauge theories, the full generaliza- 
tion of the original two-loop results for gauge groups 
with at most a single U{1) factor to the case with multi- 
ple UilVs has been formulated relatively recently, see, 
e.g., [22| and dedicated two- loop studies focusing on 
such effects in the context of, e.g., grand unified theo- 
ries (GUTs) arc available [23[. However, for the softly- 
broken SUSY gauge theories, the general two-loop evolu- 
tion equations for the soft-breaking parameters in pres- 
ence of the C/(l)-mixing effects have not yet been given. ^ 

In this study, we aim to fill this gap by presenting a 
set of substitution rules which generalize the results of 
lot 111 I to the case where the gauge group involves more 
than a single abelian gauge factor. 

The practical applications of these results are mani- 
fold. For instance, in SUSY GUTs featuring an extended 
intermediate C/(1)_r x U{1)b-l stage, see e.g. [25|, the 



C/(l)-mixing effects can shift the effective MSSM bino 
soft mass by several per cent with respect to the naive 
estimate where such effects are neglected. In principle, 
this can have non-negligible effects for the the low-energy 
phenomenology. In this respect, let us just mention that 
the theories with a gauged U{1)b-l surviving down to 
the proximity of the soft SUSY-breaking scale have be- 
come rather popular recently due to their interesting im- 
plications for the R-parity and the mechanism of its spon- 
taneous violation 



26l428l|. for Leptogenesis 



etc. 



This work is organized as follows: In Sect.|TI]we reca- 
pitulate the salient features of gauge theories with sev- 
eral abelian gauge factor focusing namely on the different 
renormalization conventions. A specific scheme in which 
the desired generalization of jlO| can be carried out in 
a particularly efficient way is identified. In Sect. IIIII the 
relevant substitution rules upgrading those in [10| to the 
most general form are given and the methods for resolv- 
ing some ambiguities emerging throughout their deriva- 
tion are briefly commented upon. In Sect. HVl we discuss 
illustrate the the importance of the kinetic mixing effects 
in a pair of specific models, focusing namely on the com- 
parison between the "rotated basis" method advocated 



in [10| and the full-fledged two-loop treatment. Then we 
conclude. For the sake of completeness, we add a set 
of appendices: some technical details of the renormaliza- 
tion scheme definition are given in Appendix |X1 the basic 
formulae of [lO| for a simple gauge group and their gen- 
eralization to the case with product groups can be found 
in Appendix [Bj In Appendix \C\ the interested reader 
can find details of the derivation of our main results pre- 
sented in Sect. IIIII Finally, Appendix [D] is devoted to 
several remarks on the gauge and gaugino matching in 
theories with multiple U{1) gauge factors. 



II. METHODS 

As mentioned in the introduction, going from a single- 
U{1) to the multiple-C/(l) case is not straightforward as 
it generally amounts to a qualitative change in the La- 
grangian. In particular, there is a need for an extra set of 
counterterms which, in simple words, keep the renormal- 
ized off-diagonal two-point Green's functions in the gauge 
sector finite. This, however, implies that the renormal- 
ized Lagrangian must contain a structure connecting the 
field tensors associated to different [/(l)'s in the gauge- 
kinetic terms, namely 



f^kin. 3 4-^/11/S-^ 



(1) 



where the different field tensors have been grouped into 
an n-dimensional vector F^,^ (with n denoting the num- 
ber of independent gauged t/(l) factors) and ^ is an n x n 
real and symmetric matrix. This amounts to ■^''^{n ~ 1) 
extra dynamical parameters. These quantities are then 
governed by a new set of evolution equations which have 
to be added to those governing the individual gauge cou- 
plings and other relevant parameters such as Yukawas 
etc. This, indeed, is the method adopted in some of the 
first studies of the subject, see, e.g., [22|. 

Alternatively, one can work in a renormalization 
scheme in which the ^-term in Eq. ([T]) is transformed 
out by a suitable redefinition of the gauge fields, namely. 



A -> S}''^A 



(2) 



^ In ref. uM the effect of the mixings of several U{l)s has been 
taken into account in the anomalous dimensions of the supcrficlds 
and in the beta-functions of the gauge couplings which serve as 
basis for the corresponding parts in the RGEs of the soft SUSY 
breaking parameters. 



which also leads to the canonical normalization of the 
gauge fields. This, indeed, affects the interaction part of 
the covariant derivative 



QjGA^QjGC^/^A, 



(3) 



where G is the original diagonal matrix^ of n individual 
gauge couplings associated to the n abelian gauge fac- 
tors and Qi is the vector^ of the relevant C/(l) charges. 
Similarly, the gauge-kinetic counterterm is transformed 



zT^bzT i ^ c"''zriBzri 



l/2zl/2^„zl/2;r-l/2_i =8Z^, 



(4) 

1/2 . 

where the subscript B denotes bare quantities and Z J^ is 
the original (diagonal) gauge-field renormalization factor 
Ab = Z J^ A. Hence, the ^~^'^ factor can be subsumed 
into a new set of ^ji{n — 1) "effective" gauge couplings 
whose combinations populate the off-diagonal entries of 
an "extended gauge-coupling matrix" 



G = GC^/^ . 



(5) 



and a suitably redefined gauge-kinetic counterterm. 

Thus, in this scheme, the off-diagonality in the gauge- 
kinetic part of the renormalized Lagrangian is absorbed 
by the covariant derivative, while the gauge-kinetic coun- 
terterm SZ^ is naturally off-diagonal in order to absorb 
the divergences in the off-diagonal two-point functions. 
Moreover, the simple QED-like relation between the bare 
and renormalized abelian gauge coupling matrices (omit- 
ting the tildes) 



Gr = GZ 



-1/2 



(6) 



remains intact because the relevant Ward identities that 
lead to the cancellation of Z^ and the Zq factors, c.f., 
Eq. (|A9[) . follow from the gauge invariance. Therefore, 
it is sufficient to work with a matrix-like gauge-coupling 
structure forgetting entirely about the ^-origin of its off- 
diagonal entries. 

This strategy, which is entirely equivalent to the for- 
mer one with a dynamical ^, is much more suitable for 
our task because it essentially amounts to replacing all 
the polynomials including individual gauge couplings in 
[10| by the relevant matrix structures, with no need*^ to 
deal with the evolution equations for the ^ matrix not 
discussed here. 



^ with indices in the group and gauge-field spaces, respectively 
^ with a lower index assigning the corresponding matter-field 
* Obviously, no information is lost so one can obtain the relevant 
RGEs for 5 components from the ones with the matrix-like gauge 
couplings. Indeed, the number of the off-diagonal entries in ^ is 
the same like the number of independent physical parameters 
governing the off-diagonal entries of G; here one has to take 
into account the freedom to bring G into a triangular form by a 
suitable redefinition of the U{1) charges. 



This, however, is not entirely straightforward in prac- 
tice. Indeed, the commutativity of c-numbers has been 
widely used in [lO[ in order to cast their results in a com- 
pact form. Thus, one has to be very careful to avoid am- 
biguities stemming from the generic non-commutativity 
of the matrix-like G"s. Furthermore, also the abelian 
gaugino soft masses have to be arranged into a matrix 
structure Af , which brings in an extra complication. 

In doing so, an invaluable key is provided by some of 
the residual reparametrization symmetries of the renor- 
malized Lagrangian. In particular. 



Q^ 

G 

A 



OlQ^, 

OiGOl , 
O2A, 



(7) 
(8) 
(9) 



where Oi and O2 are arbitrary orthogonal matrices act- 
ing in the group and gauge-field spaces, respectively, 
leave the interaction part of the covariant derivative 
Q[GA invariant. Under the same set of transformations, 
the gaugino mass matrix is rotated to 



M 



O2MO2 ■ 



(10) 



Naturally, these symmetries must be refiected at the 
RGB level. 

Thus, for instance, only those combinations C of G 
and 7 ex ^^ QiQf that transform as C — >■ OiCO^ are 
allowed to enter the right-hand side of the renormaliza- 
tion group equation for G. However, at one-loop level, 
there is only one structure involving a third power of G 
and one power of 7 that can come up from a matter-field 
loop in the gauge propagator, namely GG'^jG, so one 
immediately concludes that 



/3"°°P oc GG^jG . 



(11) 



The proportionality coefficient is trivially obtained by 
matching this to the single-C/(l) case. This also illus- 
trates that it is more convenient to work in the scheme 
with off-diagonal G than in the scheme with a non-trivial 
^, simply because the transformation properties of G 
(which is a general real matrix) arc more restrictive than 
the transformation properties of ^ (which is symmetric) . 
However, at two loop-level this becomes more com- 
plicated because then, for instance, all gauge couplings 
including those corresponding to the semi-simple part of 
the total gauge group mix among each other and/or with 
the relevant gaugino masses. Next, different Feynman- 
graph topologies can be subsumed under the same spe- 
cific term in [id . Ill | and, hence, ambiguities must be 



resolved, which often require some amount of a "reverse 
engineering" . 

Nevertheless, as we shall demonstrate in the next sec- 
tion, all such ambiguities, if properly traced back to the 
original diagrams, can be sorted out and a clear and ele- 
gant picture emerges. 

III. RESULTS 

In this section, we shall describe the generic method of 
constructing the fully general two-loop RGEs for softly- 
broken supersymmetric gauge theories out of the results 
of 



10 



11 



relevant to the case of at most a single abelian 
gauge-group factor. For the sake of completeness, the 
relevant formulae for the cases of (i) a simple gauge group 
and (ii) the product of several simple factors with at most 
a single t/(l) are reiterated in Appendices IB II and IB 2| 
respectively. The computation has been done using the 
DR scheme defined in 



13| 



A. Notation and conventions 

The gauge group is taken to be Ga^Gb^- . .(8)t/(l)", 
where the Gjc's are simple groups. We shall use up- 
percase indices for simple group-factors only; lowercase 
indices are used either for all groups or, in some spe- 
cific cases, for U{l)s only^. As mentioned before, the 
C/(l) sector should be treated as a whole and described 
in terms of a general real nx n gauge-coupling matrix G, 
an 71 X n symmetric soft-SUSY breaking gaugino mass- 
matrix M and a column vector of charges Qi for each chi- 
ral supermultiplet $,;. Notice, however, that Vi = G^Qi 
for each i are the only combinations of Qi and G which 
appear in the Lagrangian and, thus, all the general RGEs 
can be, in principle, written in terms of V^s and M only. 
We shall follow this convention with a single exception of 
the evolution equations for the gauge couplings which are 
traditionally written in terms of dG/dlogt rather than 
dF/dlogi - indeed, in this case we shall adhere to the 
usual practice. As a consequence, we expect an isolated 
G popping up in these equations. 



^ This will be evident from the context; we follow as closely as 
possible [Tot and when quoting results contained therein, the a 
and b indices go over all groups (simple and f/(l) groups). On 
other occasions, when referring to particular components of the 
(7(l)-related G, M and V matrices and vectors, a and b stretch 
over the U{1) groups only. 



Before proceeding any further we shall define some of 
the expressions that are used in the RGEs: 

• Ca (i): Quadratic Casimir invariant of the repre- 
sentation of supcrficld $i under the group Ga', 

• G (Gq): Quadratic Casimir invariant of the adjoint 
representation of group Ga ; 

• Sa (i): Dynkin index of the representation of super- 
field $i under the group Ga', 

• da (i)'. Dimension of the representation of $; under 
the group Ga] 

• d{Ga)'- Dimension of group Ga', 

• Sa (R)'- Dynkin index of group Ga summed over all 



chiral supermultiplets - Sa {R) = X^s 



Sai£l. 



Sa {R) Gb {R): Defined as J2 
SaiR)V^VR: Defined as ^,. 



'■ da(i) ' 
~~dJT) ' 



. Sa (R) V^MVr: Defined as E, ^'^'iSr^"' ; 

In addition, sometimes one has to deal with the explicit 
representation matrices of the gauge groups (denoted in 
jlO| by tj -'). Notice that here A is not a group index but 
rather a coordinate in the adjoint representation of the 
corresponding Lie algebra, (e.g., A = 1,..,3 in SU{2), 
A= l,..,8in5C/(3) etc.). 

Naturally, whenever we refer to results of refs. 10, 111 | 
for a simple gauge group (collected in Appendix IB ip . 
the a and b indices will be omitted. In all cases, repeated 
indices are not implicitly summed over. 

B. Constructing the general substitution rules 

Let us now sketch in more detail the general strategy 
for up grad ing the "product" substitution rules of Sect. Ill 
in ref. [10[ to the most general case of an arbitrary gauge 
group. For sake of simplicity, we shall focus on a limited 
number of terms here; the interested reader can find a 
more elaborate exemplification of the basic procedure in 
Appendix [Cl 

Let us begin with, e.g., the term g^G (r) appearing for 
instance in Eq. IB3I and, subsequently, in the substitution 
rules of [10| for product groups, Eq. IB34I It is clear that 
this has to be replaced by ^^5^ G^ (r)-|-'t/(l) part'. For 
a single f/(l), g'^G (r) = g'^y'l - KK- so this '[/(I) part' 



can only take the form^: Vj'Vr = Q'^GG'^Qr- There 
is no other way to obtain a number from two vectors 
Vr- Remarkably, this expression sums automatically the 
contributions of all the [/(l)'s. 

Similarly, Mg^G (r) (in Eq. IB15I for example) is re- 
placed by J2a Mao^Ga (r) + 't/(l) part'; the ingredients 
for the construction of the 't/(l) part' are two vectors Vr 
and the gaugino mass matrix M. Only V^MVr forms a 
number. 

In fact, this simple procedure allows us to general- 
ize many of the terms in the RGEs of 10|, lUl, Sect. II 
(and/or Appendix IB 1[) . As a more involved example, 
consider for instance the g'^tj -'Tr [t^G (r) m^] structure 
popping up in Eq. IB25I It is not difhcult to see that all 
terms where the representation matrices t appear ex- 
plicitly are only relevant for the abelian groups. Hence, 
if for a single C/(l) one has gHf^Ti [t^G (r) m^] = 
g'^S^Hi^ Up (^m^) , it can be immediately deduced 



-> 



that, in the general case, gH ■ ■'Tr [f^G {r)m'^'\ 
SiJ:AV^Vp)iVp^Vp){m%. 

The RGEs of G and M represent a bigger challenge, 
because they are matrix equations (i.e., the gauge in- 
dices remain open). On the other hand, this should be 
viewed as an advantage because all the relevant equations 
must then respect the reparametrization symmetries ([7])- 
pH)) . In this respect, let us reiterate Eqs. ©-([S]) which 
imply that V^s transform as T^ — > 02^i- These symme- 
tries are especially powerful in the /3-functions for the 
gauge couplings which, due to Eq. (|8]), inevitably take 
the generic form GVi (■ ■ ■ ) V^ for some chiral indices i,j. 
For example, g^S (i?) ^ g^ ^ j/p can only take the form 

Concerning the gaugino soft masses M, let us for in- 
stance take a look at the 2g^S (R) M term appearing in 
Eq. (|B2ip . Its generalized variant should be, obviously, 
built out of a pair of Vp vectors and the M matrix. How- 
ever, there are only two combinations of these objects 
that transform correctly under O2, namely, MVpVJ and 
VpVl M . Thus, due to the symmetry of M . one obtains 
2g^S {R) M ^ M Y^p VpV^ + E^ V^V^M. 

Another important ingredient of the analysis is pro- 
vided by the existing substitution rules linking the case 
of a simple gauge group (Sect. II in [lO[ and/or Ap- 
pendix lBT]) to the settings with group products (Sect. Ill 



in [lO| and/or Appendix IB2p . Consider, for example, 
the of S" (i?) G (R) term in (Eq. IB6P which, according 
to [101, gets replaced by J2b9l9b^aiR)GbiR), sec for- 
mula (|B3ip for the product groups. Let us recall that 
the expression 5* (i?) G (R) has a very particular mean- 
ing - it is the sum of the Dynkin indices weighted by the 
quadratic Casimir invariant, so J2b9a9b^a (R) Gb (R) = 
J2b p 9a9b " d (p) ■ With this in mind, whenever a refers 
to the abelian part of the gauge group, one should replace 
glSaip) ^ GVpVp^, Y.b9lCb[p) ^ Es5|Cb(p) + 
VpVp and da (p) = 1. Therefore, for the abelian sector, 
g'S {R) G (R) ^ Ep GVpVp^ [Eb 91Cb [p) + Vp^Vp] . 

However, sometimes even a detailed inspection of the 
underlying expressions docs not admit for an unambigu- 
ous identification of its generalized form. Then, a careful 
analysis of the structure of the contributing Feynman di- 
agrams is necessary. Remarkably, the number of such 
singular cases is rather limited and can be carried out 
rather efficiently, as shown in Appendix [Cl 



C. List of substitution rules 

Depending on the group sector (abelian or simple) , we 
get different RGEs for the gauge couplings and the gaug- 
ino masses. The parameters are then either the matrices 
G, M or the numbers gA, Ma- For the abelian sector, 
one obtains: 



C7(G)^0, 

g'SiR)^GY,VpV;f 

p 

's (R) GiR)^J2 ^^^^v [E ^^bCb (p) + y^y. 



(12) 
(13) 



g^c(fc) 

d{G) 



Gy,y^ 



k^k 



(14) 
(15) 



2g'^S (R) M -> M ^ VpV^ + ^ VpV^M (16) 



g^G{k)^VkV^ 



2g'^G (k) M ^ MVkV,[ + VkV^ M 



(17) 
(18) 



le/S* (R) C (R) Af ^ ^ i 4 {MVpV^ + VpV^M) 



Y.9lCB{p) + V;fVp 



+8Vpv; 



J2MB9lCBip) + V;fMVp 



L B 



For a simple group factor Ga, the substitution rules of [lOJ do not need to be changed except for two 



cases: 



(19) 



g'S{R)C{R)~^glSA{R) 



J2glCB{R) + V^VR 



Wg^S (R) C (R) M -^ SMaSa (R) 



J29lCBiR) + V^VR 



L B 



■ SSa (R) 



Y,MBglCB{R) + ViMVR 



(20) 



(21) 



As for the rest of the parameters in a SUSY model, the relevant substitution rules read: 



g^C{r)-^Y.9ACA{r) + V^Vr 

A 

Mg^C (r) ^ Y. ^^aqICa (r) + V^ MVr 

A 

M*g^C (r) -> ^ MXg\CA (r) + F^M+K 

A 

MM*g^C (r) ^ ^ MAM^g^CA {r) + V^MM'^Vr 

A 

g^C{r)S{R)^Y.9\CA{r)SA{R)+Y.{K''Vpy 

A p 

Mg^C (r) S{R)^Y1 ^^Ag\CA (r) Sa (R) + ^ {K^ MVp) (V^^Vp) 

A p 

g'C' (r) ^ ^ g^g|C^ (r) Cb (r) + 2 J] 5!^^ (0 (V^Vr) + (K^K)' (28) 

A,B A 

Mg^C^ (r) ^ ^ MAg\glCA [r) Cb [r) +Y, g\C A{r)[M a{V^ Vr) + {V^ MVr)] + (K^A/V;)(Vf K-) (29) 

A,B A 

g^C{G)C{r)^Y.9AC{GA)CA{r) 

A 

Mg^G (G) C{r)^Y. ^'^Ag\G {Ga) Ga {r) 

A 

MM*g^G (G) G{r)^Y. MAM\g\G {Ga) Ga {r) 



(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



gHf^Tr{t^.^^)^5lY.{V^Vp){mX 
p 

gHt^{t^m%^Si{V,^V;)imX 

ghf^Tr [t^G {r) m^] ^ Sf ^ (V^'Vp) {V^Vp) {m^ 

p 

g'G{i)T,[S{r)m^]^Y.9ACA{i)T:ASA{r)m^] +E(^:'^^p)'M 

A o 



(30) 
(31) 
(32) 
(33) 



(34) 
(35) 

(36) 



24/MM*C (i) S (i?) ^ 24 ^ giM^AflC^ (i) Sa (i?) + 8 ^ [{V.^MVp) {K^MWp) 



+ {V,^MM'<Vp) {V,^Vp) + (V^MHlVp) {V,^Vp)] 
48g^MM*C {rf -^ ^ 9WbCa (r) Cb (r) [32MaMX + 8MaM^ + 8MbM} 



(37) 



A,B 



^9aCa (r) [32MaMX {V^Vr) + 16Ma {V^AlWr) + 16Af^ {V^^MVr) + 32 {V^MM''Vr)] 

A 

[32 (V^MM^'Vr) (V^Vr) + 16 (V^^MVr) (y^AfV^)] (38) 

I 



IV. COMPARISON WITH OTHER METHODS 
TO INCLUDE [/(l)-MIXING 

A. General discussion 

So far, several approaches to the SUSY C/(l)-mixing 
conundrum have been proposed in the literature. Let 
us take a brief look at some of them and comment on 
their limitations as compared to the complete two-loop 
treatment advocated in this work. 

As we have already mentioned in the Introduction, one 
can attempt to choose a convenient pair of bases in the 
t/(l)-charge and gauge- field spaces for which the situa- 



tion might simplify [lOl [20[. For instance, it is always 



possible to diagonalize the one-loop anomalous dimen- 
sions 



7 



= Eq^Q^ 



(39) 



by means of a suitable Oi rotation Qi -^ OiQi = Q[, see 
d?]), so that 7' = Oi"fOl is diagonal. This, of course, in- 
flicts a change on the gauge-coupling matrix G — ?> OiG. 
However, if all the relevant U{1) gauge couplings happen 
to emanate from a single point, i.e., G ex 1 at some scale, 
Oi can be passed through G and absorbed by a suitable 
redefinition of the gauge fields dU where now O2 = Oi . 
This way, the one-loop evolution of G is driven by a di- 
agonal 7' and the initial condition G oc 1 remains intact. 
Thus, no ofF-diagonalities emerge in this case and it is 
consistent to work with the usual RGBs for individual 
gauge couplings, one per each U{1) factor. 

This approach, however, is generally limited to the evo- 
lution with a complete C/(l) unification. This is very of- 
ten not the case in practice, in particular in the GUTs 
in which the hypercharge is a non-trivial linear combina- 
tion of the relevant Cartans, such as in left-right models 
based on the SU{2)l <8) SU{2)r ® U{1)b~l gauge group, 
see Sect. IIVBI Moreover, not only gauge couplings but 
also the U{1) gaugino soft masses should coincide at the 



unification scale otherwise the method fails in the soft 
sector already at the one-loop level. The point is that 
only then the generalized one-loop correlation between 
the gauge couplings and the gaugino masses 



GM-^G'^ = const. 



(40) 



ensures the gaugino mass diagonality along the unifica- 
tion trajectory. 

At the two-loop level more complicated structures such 
as higher powers of charges, gauge couplings, Yukawas, 
etc., enter the anomalous dimensions and, in general, 
there is no way to diagonalize simultaneously all the evo- 
lution equations. Though there is still a trick one can 
implement in the gauge sector if the U{1) couplings do 
not unify [20| , there is no general way out in the super- 
symmetric case for the gauginos as also discussed in [3l[ . 
Thus, a fuU-fiedged two-loop approach as presented in 
this work is mandatory and, in fact, it turns out to be 
even technically indispensable if there happen to be more 
than two abelian gauge groups as, for instance, in 
32l | and many string-inspired constructions. 
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B. Simple illustrations 

Let us illustrate the importance of the kinetic mixing 
effects in a couple of simple scenarios which exhibits all 
the salient features discussed above. 



1. One-loop effects 

a. Gauge couplings: We shall consider the one-loop 
evolution of the gauge couplings in the SUSY 50(10) 
model of ref. 



25[ in which the unified gauge symmetry 



is broken down to the MSSM in three steps, namely, 
S'O(IO) -> SUCi)c ® SU{2)l ® SU{2)r ® U{1)b-l ^ 
S'C/(3)c (8) S'C/(2)l ® C/(l)fl ® C/(1)b_l -> MSSM; the cor- 
responding breaking scales shall be denoted by Mg, Mr 



and Mbl, respectively. Further details including the field 
contents at each of the symmetry breaking stages can be 
found in ref. 



25|. 



For our purposes, it is crucial that in this model the 
ratio Mr/Mbl can be as large as 10^° and, hence, the 
C/(l) mixing effects become important. Note that even 
a short SU{3)c ® SU{2)l «> SU{2)r (g) U{1)b-l stage is 
sufficient to split the gn and the gs-L gauge couplings 
such that the extended gauge-coupling matrix G at the 
Mr scale is rather far from being proportional to the unit 
matrix. Thus, there is no way to choose the Oi and O2 
rotation matrices such that both G and 7 



/ 15/2 -1 \ 
\ -1 18 / 



(41) 



are simultaneously diagonalized. Here N — 

diag(l, -^/S/S) ensures the canonical normalization 
of the B — L charge within the 5*0(10) framework. 
Therefore, the one-loop evolution equation relevant to 
the U{\)ii (E) U{1)b-l stage has to be matrix-like and 
reads in the abelian sector 



dt 



-A- 



-1: 



(42) 



where A-^ = 47r(GG'^)-i and t = log(/i//io)/27r. 

The reason that the U{1)r <E) U{1)b~l stage can be so 
long has to do with the fact that this gauge symmetry 
is broken by neutral components of an SU{2)r doublet 
pair, namely, (1, 1, +\,-l) © (1, 1, -i, +1) = X$j © X?e 
which are full SM singlets, and as such they do not affect 
the low-energy value of ay ■ Indeed, the would-be change 
inflicted on ay^ by the presence or absence of x% © X% 
is given by 

Aayi = pI AA-1(Mbl)pf cx p|: A7py = (43) 



where p^ = (-^/S/S, ■\/2/5) are the coordinates of the 
MSSM hypercharge in the U {!) r®U {1) b- l algebra and 
A7 denotes the relevant change of the 7 matrix. There- 
fore, at the one-loop level, the position of the Mbl scale 
is not constrained by the low-energy data and, hence, 
barring other phenomenological constraints, it can be 
pushed as close to the MSSM scale Ms as desired. 

However, this simple argument works only if the U{1)- 
mixing effects are properly taken into account. Remark- 
ably, if they are simply neglected, A7 receives only diag- 
onal entries and ay (Mz) becomes a function of Mbl- 
Moreover, stretching the Mbl-Mr range to maximum, 
the erroneous shift inflicted on ay^lAIz) can become as 




FIG. 1: One-loop gauge-coupling evolution in the MRV 
model [25(. The position of the GUT scale, the unified 
gauge coupling and the intermediate symmetry-breaking scale 
Mh were chosen in such a way to fit the electroweak data 
with aY^{Mz) ~ 59.73. The close-to-zero red line in the 
Mbl — Mr domain depicts the evolution of the off-diagonal 
entries of the A~^ = 47r(GG"^)~^ matrix which, at the one- 
loop level, scales linearly with log fi. The "optical discontinu- 
ity" in Qy^ at the Mbl scale owns to the generalized matching 
condition (ID2I). 



large as 4 per-cent as can be seen by comparing figures [T] 
and[2j Alternatively, in order to retain the desired value 
of aY^{Mz), one would have to re-adjust A/r by several 
orders of magnitude, c.f., FIG. [H This, however, could 
have a large impact on, e.g., the MSSM soft spectrum 
33|, and, in more general constructions, also on Mq and 
ac, with ramifications for rf = 6 proton decay etc. 

Finally, let us note that the "rotated-basis" method 
discussed in brief in Sect. II V Al ls only partially successful 
because the gn and gB-L gauge couplings do not coincide 
at the A/r. Indeed, the value of aY^{Mz) obtained this 
way, namely, aY^{Mz) = 60.93, is closer to the correct 
value than that received with no mixing at all, but still 
some 2% off the correct value. 

b. Gaugino masses: In order to fully appreciate the 
method advocated in this work, we should look at the 
interplay between the gauge and the soft sector. For ex- 
ample, at one loop-level, a simple illustration is by equa- 
tion (|40|) which ties the gauge couplings G together with 
the gaugino soft masses M. Consequently, the bino mass 
obeys at the scale Mg 

My (Ms) = ^^^^^plm„2PY . (44) 

where mi/2 is the GUT-scale gaugino soft mass 
matrix. From equation (|44p we see that the ratio 
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FIG. 2: The same like in FIG. [T] but without the kinetic mix- 
ing efTects taken into account. With the GUT-scale bound- 
ary condition and Mr as above, the low-energy value of Qy , 
namely, ay (Mz) = 62.51 (black solid lines) differs from the 
one obtained in the full-fledged calculation by as much as 4 
percent. Alternatively, if one attempts to obtain the right 
value of aY^{Mz) by adjusting the St/(2)_R-breaking scale, 
the new M'n scale must be shifted with respect to the cor- 
rect Mr by as much as 4 orders of magnitude (in blue dashed 
lines). 

Afy(Afs)/ay(Ms) depends on whether one includes the 
mixing effects or not as already noticed in ref. [3j| . Note 
that with non-universal initial conditions, i.e. mi/2 not 
being proportional to the unit matrix, the Py1t^i/2Py 
term mixes up all entries of mi/2- Moreover, in the spe- 
cial case that the abelian gauge couplings unify, even the 
one-loop gaugino sector evolution can be fully accounted 
for by the "rotated-basis" trick. 



2. Two-loop effects 

At two-loop level our method becomes already impor- 
tant in cases with gauge coupling unification at a cer- 
tain scale. We illustrate this by taking as an example 
the model presented in ref. 27[ where an intermediate 
SU{2>)c ® SU{2)l ® U{1)y ® U{1)b-l gauge symmetry 
is assumed to originate from a grand-unified framework. 
We assume two cases: (i) full gauge coupling unification 
at 2 X lO^^GeV and (ii) a small difference of 5% between 
the two C/(l) couplings caused by possible GUT-scale 
threshold effects. In the gaugino sector we assume uni- 
versal boundary conditions in both cases, but the effect 
gets even stronger if one considers in addition threshold 
effects in the gaugino sector as well. 

The results are given in table HI Remarkably, besides 



the expected equivalence of the "rotated-basis" method 
and the full-fledged calculation at the one-loop level, the 
relevant effective hypercharge gauge coupling turns out 
to be identical to the one obtained even at two loop-level 
if exact gauge couping unification is assumed. The reason 
is, that all additional states not present in the MSSM are 
charged only with respect to U{1)b-l but are neutral 
under the MSSM gauge group. In the gaugino sector 
the first deviations show up already in this case which 
however are only at the per-mile level. In case that one 
includes also threshold corrections at the GUT-scale the 
effects are at the percent level leading to shifts in the 
masses potentially measurable already at the LHC. 

Last but not least we remark, that the effects would be 
even larger if the U{1)y would result from the breaking of 
U{1)r(^U{1)b-l as discussed in the previous example. 

V. CONCLUSIONS AND OUTLOOK 

In this work, we have discussed the structure of the 
renormalization group equations in softly-broken super- 
symmetric models with more than a single abelian gauge 
group. Indeed, with multiple U{1) gauge factors at play, 
the effects of kinetic mixing among the abelian gauge 
fields must be taken into account in order to keep the 
theory renormalizable. 

Though, formally, the evolution equations available in 
the literature do not exhibit any obvious pathologies if 
such subtleties are not taken into account, the calcula- 
tions based on these formulas are in general incomplete 
and, thus, the results are internally inconsistent. This 
is even more pronounced in the context of SUSY models 
because it affects also the evolution of the soft SUSY pa- 
rameters, in particular the evolution of the gaugino mass 
parameters. 

Remarkably enough, the issue of the U{1) mixing in the 
softly-broken supersymmetric gauge theories has never 
been addressed in full generality, even at one loop. The 
main aim of the current study was to fill this gap and 
provide a fully self-consistent method for dealing with the 
renormalization group evolution of the gauge couplings 
and the soft SUSY-breaking parameters up to the two- 
loop level. 

To this end, we have studied in detail the existing two- 
loop renormalization group equations valid for the case 
of at most a single abelian gauge factor at play given in 
ref. 10|, lul and extended these to account for the most 
general ease of a gauge group with any number of C/(l) 
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One-loop results 


Two-loop results 




No kinetic Rotated basis Complete 
mixing method RGEs 


No kinetic Complete 
mixing RGEs 


No kinetic Complete 
mixing RGEs 


gvY 

QBLBL 
gBLY , gVBL 


0.4511 0.4700 0.4700 

0.4083 0.4243 0.4243 

0. -0.0723 -0.0723 


0.4487 0.4677 

0.4070 0.4231 

0.0 -0.0725 


0.4487 0.4686 

0.4131 0.4298 

0.0 -0.0725 


9Y 


0.4511 0.4511 0.4511 


0.4487 0.4487 


0.4487 0.4500 


Myy [GeV] 

Mblbl [GeV] 

Mbly,Mybl [GeV] 


196.34 218.13 218.13 

160.83 178.67 178.67 

0.0 - 62.39 - 62.39 


185.82 207.96 

154.88 173.19 

0.0 -63.10 


185.80 208.71 

144.26 161.97 

0.0 -62.15 


My [GeV] 


196.34 196.34 196.34 


185.82 185.96 


185.80 187.04 




Exact unification 


„GUT ■, fir „GUT 

gsL =1.05 5y 



TABLE I: Low energy values of the entries of the gauge coupling and gaugino mass matrices (qab, Mab) and the properly 
fitted MSSM parameters {qy, My), c.f., Eqs. dDTJI and (fDT3ll . We have fixed the GUT scale at 2 x 10^" with ga = 0.72 and 
imposed an mSUGRA boundary condition taking mi/2 = 500 GeV. All gaugino mass parameters are in GeV. At the one-loop 
level, we compare the case with no kinetic mixing effects included, the "rotated basis" and the full-fledged calculation. At the 
two-loop level, we include the case where qy and qbl are split at the GUT scale due to threshold corrections. 



factors. 

In particular, we have argued that all the U{1) mixing 
effects can be consistently included if the gauge couplings 
and the soft SUSY-breaking gaugino masses associated 
to the individual abelian gauge-group factors are gener- 
alized to matrices and these are then substituted into 
the formulae in |ld . Illl j in a specific manner. This, how- 



ever, is a highly non- trivial enterprise, mainly due to the 
non-commutativity of the relevant matrix-like structures, 
and a number of ambiguities had to be resolved. In this 
respect, the residual reparametrization invariance of the 
covariant derivative associated to the redefinition of the 
abelian gauge fields turned out to be a very useful tool, 
yet in many cases one had to resort to a detailed analysis 
of the relevant Feynman diagrams. 

The general method has been illustrated for two cases: 
(i) at one-loop level where due to an breaking of the orig- 
inal group two C/(l) factors emerge with different gauge 
couplings and (ii) at the two-loop-level in a model where 
gauge coupling unifications occurs but where threshold 
corrections are taken into account. In both case we ob- 
tain effects in the percent range and we remark, that 
none of the previously proposed partial treatments can 
account for the full effects. 

Last but not least, let us stress again that our results 
are completely generic and, as such, they do not require 
any specific assumptions about the charges of the chi- 
ral multiplets in the theory and/or the boundary con- 



ditions applied to the relevant gauge couplings. This 
makes the framework very suitable for implementation 
into computer algebraic codes calculating two-loop renor- 
malization group equations in softly-broken supersym- 
metric gauge theories such as SARAH I35l - l37l | and Susyno 

ii. 
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Appendix A: Renormalization of QED ® QED 

In this appendix, we comment in more detail on renor- 
mahzation of abehan gauge theories, focusing on the sim- 
plest non-trivial case exhibiting the effects of kinetic mix- 
ing, namely the "QED-squared" scenario featuring two 
independent abehan gauge groups f/(l) (g) U{1). 

Let us start with the basic bare Lagrangian of QED^ 
including an explicit kinetic-mixing term 






(Al) 



Here Vj a-rc the relevant matter fields (whose number 
must be in general equal to or greater than the number 
of the abclian gauge factors otherwise there is no way to 
distinguish among all stipulated U{1) factors). A stands 
for a 2-component vector (in the group space) comprising 
the gauge fields associated to different U{1) factors, G 
stands for a (so far formal, i.e, diagonal) 2x2 containing 
the relevant pair of gauge couplings and ^ is a symmetric 
and real 2x2 matrix parametrizing the gauge-kinetic 
form. 



1. Scheme A: a non-canonical gauge propagator 
and diagonal gauge couplings 

Leaving ^b in the game, one defines the renormalized 
and the counterterm Lagrangians as 



C = TpSP - mi)ipi ■ 



5C = iip^SZ^.p-ipi - iljiSZm,m,ip. 



^iQjGU, ~ \F^u^F^^- , (A2) 

(A3) 



where 



SZgG = Z^.Qi GbZ^ — Qj^ G , 



s^ = z'J'CbZ^ 



.1/2 



^: 



(A4) 
(A5) 



and 5Zm and SZ^. arc unimportant for our considera- 
tions. These countcrterms are fixed by the renormal- 
ization conditions so that they render the renormalized 
Green's functions of the theory UV-finite. For a diago- 
nal Za and for any fixed ^b, the off-diagonal entries in 
S£, cannot be matched by the right-hand side of Eq. (jASp 
unless ^ is a dynamical quantity. Remarkably, in this 



GbCb '^b 



scheme the gauge coupling can be retained in a diagonal 
form throughout the KG evolution. This is because the 
relation between the bare and the renormalized couplings 

I/O 

Gb = GZ^ can be brought into the form (trading Za 
for ^ and S£^) 

= GZ-'^^^s'Z-'^^G^ = G(e + SO-'G^ 

(A6) 
(which holds to all orders in perturbation theory) from 
where it is clear that any non-diagonal entry of the RHS 
of the evolution equation for G can be absorbed into ^. 
So, in this scheme, ^ is a dynamical quantity while G can 
be kept diagonal. 



2. Scheme B: a canonical gauge propagator and 
non-diagonal gauge couplings 

If, instead, £, is absorbed by a suitable gauge-field re- 
definition A — > S^^/^A = A into matrix for the coupling 
constants, one is left with 



C = i^M - mO^. - V.QfGM - IF^.F^" , (A7) 
SC = iipiSZ^.-P^Ji - ip.SZmi'mi'ipi - TpiQlSZ^G^^^'i 







where 




QfSZ^G = Z^.^QjGBZf-QjG, 


(A8) 


6Z^ = (Zf)-Zf-I, 


(A9) 


with ZV^ = ^-i/2^y2^V2 ^j^^^ g^g heiore, 





Gb = GZ 



-1/2 



(AlO) 



It is again clear that the non-diagonality inflicted on Z^ 
by the renormalization conditions renders the RHS of the 
gauge-coupling evolution equation non-diagonal. How- 
ever, in this scheme, ^ has been swallowed by the gauge- 
field renormalization counterterm and, as such, does not 
need to be treated as an extra dynamical quantity. In 
other words, the whole effect is accounted for by the off- 
diagonal form of the generalized gauge coupling G. 



Appendix B: Recapitulation of the t^vo-loop RGBs 
for simple groups and their products with at most 

one f/(l) 

1. Case A: Simple gauge group 

For completeness we display here the RGEs in the case 
of a simply gauge group based on 10l4l2l| . For a general 
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A'^ = 1 supersymmctric gauge theory with superpotential with 



w{<^) = L,$, + ]^^i''^^<^J + iy''^^$,$,$fc , (Bi) 



the soft SUSY-breaking scalar terms are given by 



+ c.c. 



+ [m'YjCp.^,* + -M\a\a 



(B2) 



Here we will follow [10[ and assume that repeated indices 
are summed over. Note also that lowered indices imply 
conjugation (e.g., F^fc = Y"^^^ ). In the notation defined 
in sect. IIIIl the anomalous dimensions of the chiral su- 
perfields are given by 



„(i)j 



-y.p^y^p^ - 2<5^Vc(z) , 



(B3) 



yf ^' = 9^Y,j,,Y^p'i[2C{p) - c{i)] - iK,„„r"p«yp,,r™'^^- 

+ 25lg^[C{{)S{R) + 2C{if - 3C(G)C(z)] , (B4) 



and the /3-functions for the gauge couplings are given by 



/?W =53 [S{R) - iC{G)] 



(B5) 



Pf^ =9^ {-6[C(G)]^ + 2C(G)S{R) + 4S{R)C{R)} 

- g^Y'^''Y,,kC{k)ld{G) . (B6) 

The corresponding RGEs are defined as 



1 



Jt' 



IGtt' 



r/3, 



(1) 



(167r2 






(B7) 



Here, we used t = InQ, where Q is the renormalization 
scale. The /3-functions for the superpotential parame- 
ters can be obtained by using superfield technique. The 
obtained expressions are 



/3^' 



/3^ 



yljp 



1 



2 'P ^ 



>% 



(2)fc 



IGtt^ '^ ■ (167r2)2 'P 
+ (fc <^ i) + (fc O j) , 

1 .1^,,■ 1 



?P 



/?! 



M 



LP 



IGtt 
1 



2 'p 



2 /p ^ 



(167r2)2 

1 



(2W 



r7i 



(2)« 



(B8) 

+ (i ^ *) , 

(B9) 
(BIO) 



167r2 'P ■ (167r2)2 'p 
The expressions for trilinear, soft-breaking terms are 



dt 16^2 



(1) 



ijk 



1 



(167r2) 



2\2 



(2) 



ijk 



, (Bll) 



?(i) 



ijfc 






Y'^'Yirar.h"'^'' (B12) 



(/I'J'''^ - 2MY'^'') g^C{k) + {k ^ i) + {k ^ j) , 



(2) 



ijk 



— /yy'y, v^wgy Y 

-'^ ^Lmn^ ^ r)nr ^ 



mrk 



-* ^Imn-^ ^pqr"' ^ ^ Lmn"' ^ pqr ^ 



+ (h'^'YipgYP'^'' + 2Y'^'YipqhP'^'' 

- 2MY'''Yip,YP''''')g^ [2C{p) - C{k)] 
+ {2h'^^ - MIY'^^) g^\c{k)S{R) + 2C{kf 

- 3C(G)C(fc)l +{k^i) + {k^ j) . 



(B13) 



For the bilinear soft-breaking parameters, the expressions 
read 



dt 167r2 



P'^' 



(167r2) 



2\2 



,(2) 



(B14) 



with 



/3l^^ 



'' = ^fo^'rz„„r™"^' + ir*^''rz„„6™" 



A 



+ M''>1mn/i™"^' - 2 (6'^ - 2M^i'=) g^C{i) + {i ^ j) , 

(B15) 

^3 1, 



(2 



-^ ^ Imn ^ ^ pqr ^ 



1 



^ ^ ivanh^ ^pqr"' P ^ivan'^ ^ pqr ^ 



h*' ^Iran^ ^pqr"' ^ ^ hnn'^ ^ pqr ^ 

+ 2Y'^'Yip, (6P« - fiP''M)g^C{p) + (b^'Yip^YP"^ 
+ 2M''rip,/iP« - 2fi''Yip,YP''ni)g^ [2C{p) - C{i)] 
+ (26'J' - S^i'm) / \cii)SiR) + 2C{if 

- 3C{G)C{i)] + (i o j) , (B16) 



Finally, the RGEs for the linear soft-breaking parameters 
are 



dt 167r2 



■pf 



(167r2) 



2\2 



4^) 



(B17) 
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with 



/?«■ 



-Y^'^Ypir^SP + LPY^inh^'^ 



+ M'^^info'" + 2Y'^^\m')'nki + h'^'^hk 



(B18) 



/^r 



2.g2c(0r'^'rpfe/5P - -Y'^^^Yq^tY'^'Y^uS-P 



- Ag^C{l){Y'''^M - h'^^)YpkiL^ 

- [Y''''^Yq,th''%ki + h'''''Y^,tY''%ki]LP 

- Ag'cm.^iiti'^'M - 6"V - [Yjn.hO^'Yistt^"' 

+ Y'P\m')l^jiki+Y'''^{m%^jiki) 

Y'^^^Yq^th'^ki + h''''iYq,tY'Hki 
+ h''"'h^,tY'^'m + Y'P'i{m')%,tY'''^lu 

+ Y'^'iY^,tYP'\mXm + r^'^p(TO2)^r,,ty'^Vfc; 
+ 2Y'^'iY^sM''Tty'''m + y'''''Ksth'''m ■ 

(B19) 

With these results, the hst of the /3-functions for all cou- 
phngs is complete. Now, we turn to the RGEs for the 
gaugino masses, squared masses of scalars and vacuum 
expectation values. The result for the gaugino masses is 



dt 



M 



1 









(B20) 



with 



/3iy =g^ [2S{R) - 6C(G)] M , (B21) 

/3f^) =/ {_24[C(G)]2 + 8C(G)S'(i?) + 16S'(i?)C(i?)} M 

+ 2g2 [/i^J'fc _ MY'^^] Y,jkC{k)/d{G) . (B22) 

The one- and two-loop RGEs for the scalar mass param- 
eters read 



d 



1 



167r2 



(1) 



1 



(167r2)2 



A 



(2) 



(B23) 



with 



(1) 



' 1-. -. 2^J 1 



|y.p,rP'?"(m2); + -Y=P^Y,,^{ni% 



+ 2K,p,y^P'^(m2)^ + /i,p,/i^P? - 8(5^MMtg2c(z) 



,(2) 



- tpq J 

2.g2tf^Tr[t^m2] , 
' = -\{m')%anY"'^^Yp,rY^'>^ 

)\ /I rnri ^ ^ pgn 



(B24) 






^ilrn^ V 'n PQ^ ^ilm-^ "'npq"' 

OY, Yi'-^Y Y"^P^(m'^\'^ ~ h-, hi^"-Y Y"^P1 

V, VJ"'"l'm2'l' V V''P9 — V-, /iJ'"/7 Y™P1 
J-iini-^ V Jr ^PQ ^ihn"' if'npq^ 



2\1 
I 



■ 4Y,pgY^P''MM'' 

^2.Ajf.A^2U 



g' [C{p) + C(q) - C{i)] 



- 2gX\i''^'yrYipqY^P'^ + 8/tf^Tr[t^C(r)m1 



4.Ajr 



Sfg^MM'' 



24C{i)S{R) + ASC'iif - 72C{G)C{i) 



8<5f/C(z)(Tr[5(r)m2] - C{G)MM^) 



(B25) 



The RGEs for a VEV f * are proportional to the anoma- 
lous dimension of the chiral superfield whose scalar com- 
ponent receives the VEV 



—V' = V'' 
dt 



1 



IGtt' 



'Ip 



[1)1 



1 



(167r2)2 ^P 



(B26) 



2. Product groups with at most one L'^(l) 

To generalize the formulas above to the case of a direct 
product of gaug e groups, the following substitution rules 
are needed [lO[. Note, we give these replacements here 
only for completeness and they are not sufficient in the 
case of several f7(l) gauge groups, see sec. HII] for the 
necessary extensions. 

For the 13 functions of gauge couplings and gauginos the 
rules are 



g^C{G) ^ glG{Ga) , 

g^S{R)^glSa{R). 

g^C{Gf ^ glG{Gaf , 

g''G{G)S{R) ^ glGiGa)Sa{R) 

b 



(B27) 
(B28) 
(B29) 
(B30) 
(B31) 



g'SiR)GiR) ^Y.9l9lSa{R)Cb{R) , 

6 

lQg^S{R)G{R)M -^ 8Y,glglSa{R)Gb{R){Ma + M,) , 

b 

(B32) 



g''C{k)/d{G)^glGa{k)/d{Ga). 
For all the other f3 functions, we need 
g^Gir)^J29lCa{r), 

a 

/C(r)5(i?)^^5aCa(r)5a(i?) 



(B33) 

(B34) 
(B35) 
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/C(r)C(G) ^ Y. 9lCa{r)C{Ga) , (B36) 

a 

a b 

(B37) 

a b 

X (32MaMi + 87\/,Mt + MlbMl) , 



(B38) 
gV^Tr(t^m2) ^ ^ .g^(tf HTr(t>^) , (B39) 

a 
a 

(B40) 

a b 

X Tr[t^a(r)m2] , (B41) 
g^C{i)Tr[S{r)m^] ^ ;^,9^a(^)Tr[5a(r)m2] . 

a 

(B42) 



Appendix C: Obtaining the substitution rules 

In this Appendix, we illustrate in more detail the meth- 
ods used throughout the derivation of the substitution 
rules given in section Hill 



1. The role of the Vi vectors and the A/ matrix 

As mentioned in the text, the U{1) gauge coupling ma- 
trix G and the charge vectors Qi of the chiral superfields 
$i appear always through the combination Vi = G^'^Qi. 
The only exception are the RGEs of G, where there 
should be a leading free G. For example, the ip^'' ip'^ Af^ 
vertex is proportional to V"" (component a of the vector 



V. 



'ipl 



V, 




Aa 



(CI) 



i*iAi Aa J,*ij,i Ao. Ab Mi-*nl,i\a 



AlAl, (f)"'ilj'X'' and 



,'M*rhO 



are propor- 



Similarly the vertices 
the Yukawa independent part of 0**(; 
tional to F,", V^-V^^ , V^ and V^Vj respectively. In addi- 
tion, there is to consider the [/(I) gaugino mass matrix 



M, 



Mab 

where Mab is the a, b component of M . 



(C2) 



2. RGEs with no f/(l) indices 

Only diagrams underlying the RGEs for G and M con- 
tain external U{1) gauge bosons/gauginos. As such, in 
all other equations, while vectors Vi and the matrix M 
may be present, they must be in combinations that are 
scalars with no free U{\) indices. 

Consider Mg^C (i) appearing in the one-loop RGE of 
the bilinear scalar soft terms 6*^ , which is to be replaced 
by J2a MAg\CA {i) + Vi^MVi. The simple groups con- 
tribution does not interest us though, so we shall neglect 
it. We can see that V^'^MVi is the only structure that 
can generalize the expression Mg^C (z) = Mg^yf for one 
U{1) group only. Observe also the contraction of the 
t/(l) indices in the expression - it comes from the possi- 
bility of having any of the U(l) gauginos in the internal 
lines of the contributing diagram. 



(C3) 



The amplitude is proportional to J2a b ^t^abV^ ^"^^ = 
fi^^Vi^AIVj. Note that, for any pair of values i,j the 
gauge symmetry forces ^'-^ = unless Vi + Vj — which 
means that n^^Vj = —fi^^Vi so the amplitude of the dia- 
gram is indeed proportional to VJ^MVi. 

This requirement that expressions with V^s and M's 
must form scalars is enough to derive the Eqs. (PO)) - 
(HH), dm-dMl) and dSHl) from the existing substitu- 
tion rules for gauge groups with multiple factors. We 
are left with the terms g'^C{r)SiR), Mg^C{r)S{R), 
5''C (i) Tr [5 (r) m2] and 2Ag'^MM*G {i) S {R). Note 
that one can write S [R) as Tr [S (r)] in the notation of 
ref. jlO[, so in all four cases there is a sum over field 
components of chiral superfields. For diagrams with up 
to two-loops and with no external gauginos nor gauge 
bosons, the factors S (i?) and Tr \S (r) tti^] can only come 
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from the following sub-diagrams'': 




r\/\/w 



(C4) 



rwx/v 



Take for example g^C (r) S (R) - Y^pd^lvl- The 
reason why one cannot immediately generalize this ex- 
pression to include U{1) mixing effects is because in 
theory it could take the form ^ {Yt^^p) (K"^^) oi' 
T.p{^^'^r){V^Vp). But looking at the above dia- 
grams, such ambiguities go away because in all cases 
the y's which are summed over (the V^'s) do not con- 
tract with each other. They contract with something 
else at the other ends of the gauge boson/gaugino 
lines. As such, there are no V^Vp^s, in these expres- 
sions and with this piece of information, combined with 
the known rules for a gauge group with multiple fac- 
tors, Eqs. dMl), (Ell) and ^ follow. The substitu- 
tion rule given in Eq. dST]) for 24,g^ M M* C {r) S {R) 
appearing in the two-loop equation of the soft scalar 
masses is more complicated since the placement of the 
M, M'l' gaugino mass matrices between these Vs is rel- 
evant. Nevertheless, from the following diagrams we 
can calculate that the C/(l)'s contribution to this term 
is 8Ep [{V^'MVp) {V^'M^) + {V^'MM^) {V^'Vp) 




(C5) 



^ It is conceivable that they could come also from diagrams with 
one (f)* <!>* (jxj> vertex, but we may choose an appropriate gauge, 
the Landau gauge, where these are because an external scalar 
line always couples to a gauge bosons at a three point vertex. 




(C6) 



(C7) 



3. RGEs with [/(I) indices 

The RGEs for G and M are the only ones with free 
C/(l) indices. For the beta functions of the gaugino 
masses, we will be interested in looking at diagrams with 
two incoming gauginos. As for the coupling constant, 
due to the Ward identities, the contributing diagrams 
are those with two external gauge bosons. From the am- 
plitude of these diagrams we still have to add a G factor 
in order to obtain fie- Pictorially, 



EGaa' X 




(C8) 



Note that all the terms in Pq must be of the form 
GV^ {■ ■ ■)Vj for some i,j as mentioned in the text. 
These qualitative considerations suffice for our purposes. 
Keep also in mind that 

(i) the RGEs arc invariant under the set of transfor- 
mations G -^ OxGOl; V^ ~f O2V,, M -^ O2MOI 
for any orthogonal matrices Oi, O2] 

(ii) M is a symmetric matrix and so must be ^^. 

Taken together, these considerations allow us to deduce 
Eqs. dHJ-dUl) (Eq. ^ is trivial). 

We shall exemplify this for the case of 
\&g'^S {R) G (R) M which for multiple factor groups 
is replaced by 8J2^glgiSa{R)Cb{R) {Ma + Ah) 
in the RGEs of Ma- This is the same as 

8j:p,i,9l9 f''ffpfH Ma + Mb). Groups a and b 
are independent so the expressions ^b9b^b{p), 



J2b Mbgb^b (p) are decoupled from g'^ 



2 Sa(p) 

dlipj 



Ma9i 



2 Sa(p) 
daip) ■ 
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Inclusion of C/(l) mixing effects in the first pair of 
expressions is easy because there are no free C/(l) 
indices: Y^bdlCbip) -> Y. b 9lC b {p) + Vj Vp and 
Eb Mt^glCi, {p) ^ Eb MBg%CB ip) + Vj'MVp. If the 
group a is an C/(l), then in case of a single C/(l) this 
corresponds to g^ j^l^l = g^Vp which generahzes to 

9a d''() ~^ ^p^p ■ Similarly, the only symmetric matrix 
expression which respects the O2 symmetry that can 
generalize Afaff^f^ is \{MVpV^ + VpV^M). 

Assembling these pieces gives Eqs. [19] for 
lQg^S{R)C{R)M. The structure of the final ex- 
pression is verifiable by looking at the relevant diagrams: 




^<^v^v ^^v^^v^ l 



^A^r>A^^ i 







5> y 



Appendix D: Gauge-coupling and gaugino-mass 
matching 

In this appendix we comment on yet another key ingre- 
dient of any practical application of the methods advo- 
cated in this work, namely, on the matching between the 
symmetric and asymmetric phases of spontaneously bro- 
ken models with multiple U{\) gauge groups, where the 
latter is described by an effective theory with a reduced 
dimensionality of the abclian sector. 



1. Gauge couplings 

Whenever a charged field develops a VEV, the original 
t/(l)a ® U{l)b gauge symmetry gets broken down to a 
single U{l)c spanned over the unbroken combination of 
the original generators 

Q1 = PaQ'i +PbQ\ (no summation). (Dl) 

A gauge coupling associated to this residual symmetry is 
then in general a function of the original gauge couplings 
ga and gti, the p-coefHcients above, and, at higher loop 
order, also other quantities such as group Casimirs etc. 
Let us recapitulate in brief how the matching between 
the two regimes determines gc, focusing mainly on the 
"tree-level matching" as needed when working with the 



one-loop RGEs. We comment on the changes for two- 
loop evolution at the end of this section. 

It is convenient to address the issue in two steps. First, 
one can consider the matching between two QED^ sce- 
narios which are connected just by a pair of O-rotations 
as in Eqs. d!!)-®. In particular, any specific choice of Oi 
represents a transition from the original set of the U{1) 
charges Q, = [Q'i^Q'iY ^ a new set g', = {Q^„QfV 
where Q'^ = OiQi- This, in turn, transforms the rele- 
vant (matrix of) gauge couplings as G ^- G" = OiG so 
that the interaction part of the covariant derivative re- 
mains intact, QfGA — Q'^ G'A. Regardless of whether 
simultaneously an O2 rotation has been performed on 
the gauge fields, i.e., A ^ A' = O2A (inducing G — > 
G' = OiGOl), one always has G'G''^ = dGG^Of , or, 
equivalently. 



{G'G")-^^Oi{GG')-^Oi 



(D2) 



which, as we shall see below, is more useful in practice. 
Hence, for any specific choice of Oi, this equation yields 
the link between the specific non-linear combinations of 
gaa, gab; gba and gbb entries of G and the gcc, gcd, gdc 
and gdd entries of G' and, as such, provides the desired 
matching condition. Note that, given the Pa and ph co- 
efhcicnts in Eq. (|Dip . the first row of Oi is fixed and 
its second row is determined from orthogonality up to a 
global sign. 

However, O2 plays an important role if one, e.g., needs 
to get a coupling associated to a specific gauge field, for 
instance the one associated to the "residual" U{l)c that 
survives the U{l)a^U{l)i, breakdown. This gauge boson 
corresponds to the massless eigenstate of the gauge-boson 
mass matrix 



Ml = G^{HyQ*HQji{H)G 



(D3) 



where {H)is theU{l)a<E) C/(l)6-breaking VEV. Since the 
conserved charge Q^ in Eg . ID 1 1 annihilates this VEV, it is 
convenient to go into the primed basis where Q' jj{H) = 
(0,T/). Thus, 



Mi 



G'^01{H)^Q'*hQh{H)0iG 

'q 

V^ 



TnT 



G'0{ 



OiG 



(D4) 



A convenient gauge-field transformation A ^f A' — O2A 
would bring this mass matrix into a diagonal form 



Ml 




(D5) 
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(using G" = O1GO2) if and only if G' is upper-triangular, 



I.e. 



G' 



g'cd 



g'dd 



(D6) 



The zero in the 12 entry has a clear physical implication: 
for the couplings of the surviving (massless) A'^ gauge 
boson only g^ = (?^c i^ relevant in the effective theory 
whereas one does need to care the g^^ and g'^^ couplings 
of the heavy A'^ which is integrated out except for the 
calculation of higher-dimensional effective operators. 

From here it is also easy to understand why Eq. ()D2p is 
better suited for practical purposes than the one without 
inverse: Indeed, the 1-1 entry of (G'G""^)^^ reveals g'^^ in 
a very simple way, namely. 



'r<'T\ 



{G'G- 



(5^ 



9cd9dd 9cc 

+ 9?d)9'd-d'9'c-c' 



(D7) 



(the dotted term follows from symmetry) so that it is 
sufficient to look at the 1-1 entry of the RHS of Eq. (|D2[I . 
whilst without the inverse 



G'G 



IT 



11 I I 
9cd 9cd9dd 

9'Id 



(D8) 



and in order to extract g'^^ one has to solve a non-linear 
system involving all three independent entries of the RHS 
of Eq. dm]), namely, OiGG'^Oj. 

At higher-loop orders, the situation becomes slightly 
more involved, especially if the [/(I) (g) U{1) gauge struc- 
ture is tensored with a semi-simple gauge factor Gx- 
For example gauge-boson canonical normalization effects 
have to be considered at the two- loop level 39|,|40[, which 
yields for example extra group-Casimir factors associated 
to Gx entering formulae like Eq. (|D2p . see, e.g., [23|. The 
specific shape of these terms is, however, renormalization 
scheme dependent. 



block-diagonal form, c.f., Eq. (jD5[) . and with Oi rotation 
imposed on charges) 



M 



\',H 



( M' V2G'^Q'h{H)* V2G'^Q'h{H^)*\ 
Whh Whh- 



I 

(D9) 

Here M' = O^MO^ is the gaugino soft mass matrix, 
W is the superpotential; the subscripts of W denote the 
derivatives of W with respect to the superfield H and 
its charge conjugate W^ . As before, due to the specific 
choice of the primed basis one has Q']j{H) = (0, V) and 
Q' jj{H''') = (0, y*). Furthermore, the triangular shape of 
G", see Eq. (jD6|) . ensures that the Ac gaugino correspond- 
ing to the first row/column receives no mass contribution 
due to the spontaneous symmetry breaking. 

Thus, one concludes that the effective soft mass of the 
gaugino associated to the surviving gauge group C/(l)c is 
given by the c-c entry of the M' mass matrix 



Mr = ML 



(DIO) 



Note that, in principle, there is no need to calculate the 
O2 matrix because one can trade it for the gauge cou- 
plings and the (known) Oi: 



O2MOI = G'^OiG-^^MG-^OjG' 



and thus 



G'-^'^M'G'-^ = OiG-^'^MG-^Ol 



(Dll) 



(D12) 



The main advantage of this formula is that, again, the 
LHS reveals {M')cc in a particularly simple manner, 
namely 



G' 



■/-IT 



M'G'-^ = 



MU9'c 



(D13) 



2. SUSY and the gaugino masses 

Concerning the gaugino masses, the situation is only 
slightly more involved here. Sticking to the simplest 
U{l)a ® U{l)b — > ^(l)c case as before, the gaugino(A)- 
higgsino(-ff) mass matrix reads schematically (in the O2- 
rotatcd basis bringing the gauge boson mass matrix to a 



with R12 = -9cdKc + 9'ccKd9'dd/9cc and R22 = 
i9clM:id-29cc9cdKd + 9?dKc)/92d9cl and, thus, in 
combination with (JD7|, admits for a simple extraction 
of the surviving gaugino effective soft mass. Moreover, 
due to the one-loop RGE invariance of the G^^^MG^^ 
combination, the RHS of Eq. (jD12|) is directly connected 
to the initial condition. 
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